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We analyze the evolution of two-band superfluidity from the weak coupling Bardeen-Cooper- 
Schrieffer (BCS) to the strong coupling Bose-Einstein condensation (BEC) limit. When the inter- 
band interaction is tuned from negative to positive values, a quantum phase transition occurs from 
a 0-phase to a vr-phase state, depending on the relative phase of two order parameters. Furthermore, 
population imbalances between the two bands can be created by tuning the intraband or interband 
interactions. We also find two undamped low energy collective excitations corresponding to in-phase 
and out-of-phase modes. Lastly, we derive the coupled Ginzburg-Landau equations, and show that 
they reduce to coupled Gross-Pitaevskii equations for two types of bosons in the BEC limit. 



PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 

A two-band theory of superconductivity was intro- 
duced by Suhl et al. 1] in 1959 soon after the Bardeen- 
Cooper-SchriefFer (BCS) theory to allow for the possibil- 
ity of multiple band crossings at the Fermi surface. This 
model has been applied to high- Tc superconductors and 
MgB2, where, in the latter case, experimental properties 
can be well described by a two-band weak coupling BCS 
theory |^|^^|^. Unfortunately, interband or intraband 
interactions can not be tuned in these condensed matter 
systems, and their properties can not be studied away 
from the BCS regime. However, two-band fermions may 
also be produced experimentally with ultracold atomic 
Fermi gases in optical lattices p or in single traps of 
several hyperfine states. In this case, (intraband and in- 
terband) interactions may be tuned using Feshbach reso- 
nances which allow for the study of the evolution of two- 
band superfluidity from the BCS to the Bose-Einstein 
condensation (BEC) limit. The BCS to BEC evolution 
in the two-band problem is much richer than the one- 
band case which has already been experimentally stud- 
ied IE S IS 0j IQij >1L il^j: since additional interaction 
parameters (interband and intraband) may be controlled 
externally. 

Furthermore, a two-band model may be used not only 
to describe experiments involving several hyperfine states 
of the same fermion, but also may be used to study two 
different fermionic species (e.g. ^Li and ^"K). In systems 
involving mixtures of two different alkali atoms, simple 
one-band theories may not be sufficient to describe the 
interactions between the two species of atoms, and two- 
band theories may be necessary to model future experi- 
ments. Thus, due to recent developments and advances in 
atomic physics described above, and in anticipation of fu- 
ture experiments, we describe here the BCS to BEC evo- 
lution of two-band superfluids for all coupling strengths 
at zero and finite temperatures. 

The main results of our paper are as follows. We show 
that a quantum phase transition occurs from a 0-phase 
to a TT-phase state (depending on the relative phase of the 
order parameters of the two-bands) when the interband 



interaction J is tuned from negative to positive values. 
We found that population imbalances between the two 
bands can be created by tuning intraband or interband 
interactions. In addition, we describe the evolution of two 
undamped low energy collective excitations correspond- 
ing to in-phase phonon (or Goldstone) and out-of-phase 
exciton (finite frequency) modes. Near the critical tem- 
perature, we derive the coupled Ginzburg-Landau (GL) 
equations for a two-band superfluid, and show that they 
reduce to coupled Gross-Pitaevskii (CP) equations for 
two types of composite bosons in strong coupling. 

In order to obtain the results described above, we start 
from a generalized Hamiltonian for multi-band superflu- 
ids with spin (pseudo-spin) singlet pairing 



H = ^ 6i,£7(k)at^^(k)a„,^(k) 

n,k,cr 

- E KI™6L(q)&..(q), (1) 

n,m,r,s,q 

where the indices n,m,r and s label different en- 
ergy bands (or components), and a labels spins (or 
pseudo-spins) . The operators ^ (k) and (q) — 
Ekr;™(k)ai,^(k + q/2)al j-k + q/2) create a single 
and a pair of fermions, respectively. The symmetry fac- 
tor r„m(k) characterizes the chosen angular momentum 
channel, where r„m(k) = knm.o/ {k'^m + ^^)"'^''^ is for the 
s-wave interaction in three dimensions. Here, knm,o ~ 
Rnm.o s^ts scale at small and large momenta, where 
Rnm,a pls^ys the role of the interaction range. In addition 
^„,o.(k) = e„(k) - ^„^cr, where e„(k) = e„,o + A;^/(2M„) 
is the kinetic energy {h = 1) and M„ is the band mass of 
the fermions. 

From now on, we focus on a two-band system such that 



yr 



Vnr^nm^rs with distinct intraband (Vlij^22) > 



and interband (Vl2 = V21 = J) interactions. Notice 
that, J plays the role of the Josephson interaction which 
couples the two energy bands. In addition, we assume 
that the total number of fermions is fixed {N = Ni + 
N2) such that the chemical potentials of fermions are 



2 



identical {pin,a ~ /^), and that the reference energies are 
such that £1^0 = and e2,o = e_D > 0, as shown in Fig.^. 
Here, eo = fc|,/(2Mi) < ep ^ kl p/{2Mi), where ep 
is the Fermi energy and fc„,_F are the Fermi momenta 
fci,F = kp and fc2,F — [2M2{eF — ^d)]^^^- Since the low 
energy physics depends weakly on fc„^o in dilute systems 
{AfRg <C 1) characterized by fc^ q ^ _f where M = 
N/V is the density of fermions and V is the volume, we 
assume for simplicity that fc„^o = fco S> kp. 
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FIG. 1: Schematic (a) figure of two bands with reference 
energies ei,o = and £2,0 = ec, and (b) phase diagram of 
0-phase and 7r-phase states. 

The gaussian action for H is {ks — 1,P — l/T) 

^gaus. = So + At(-q)F-i(g)A(g), (2) 

where q = (q, ive) denotes both momentum and bosonic 
Matsubara frequency = 2£tt/P. Here, the vec- 
tor At(-q) = [Al{q), A^{^q), A^i-q), A^(g) ] is 
the order parameter fluctuation field, and the matrix 
F~^{q) is the inverse fluctuation propagator. The sad- 
dle point action is Sq = -/? ff«m A* oA,„,o + 
E„,k {/3Kn(k)-£;n(k)]-21n[l + exp(-/3£;„(k))]}, where 
.911 = -V22/ detV, 522 = -Vii/ det V and 512 = 321 = 
J/detV with det V = 1/11^22 - > 0. Here, Ej,k) = 
[^^(k) + |A„(k)p] 2 is the energy of the quasiparticles and 
A„(k) — A„^orn(k) is the order parameter. 

The action given in Eq. Q leads to the thermody- 
namic potential figauss = f^o + f^fluct, where flo — Sq/ (3 is 
the saddle point and fifluct = /3"^ jZq lndet[F-i(q)/(2/3)] 



where the matrix elements are given by Onm = —gnm — 
<5«mEk|rm(k)ptanh[/3^,„(k)/2]/[2S,„(k)]. Here, (5„™ 
is the Kronecker delta. Notice that the order parame- 
ter amplitudes are the same for both the 0-phase and 
TT-phase as can be shown directly from Eq. but 
their relative phases are either or tt. In what fol- 
lows, we analyse only the 0-phase state, keeping in mind 
that analogous results (with appropriate relative phase 
changes) apply to the 7r-phase state. We can ellimi- 
nate Vnn in favor of scattering length a„„ via the re- 



lation 1/Vn- 



-M„V/(47ra„„) + |rn(k)| V[2e«(k)], 



which can be solved to obtain 1/ {kn^pann) — ko/kn^p — 

^TT/ikn^pVnnM^V). 

The order parameter equation needs to be solved self- 
consistently with the number equation N = —dQ/dj-i 



leading to 



gauss 



N, 



Nq + A^fluct, and is given by 

a[dctF-i(q)]/9/i 



gauss 



detF-i((j) 



(4) 



Here, the first term is the saddle point (A^i -I- A^2) and 
the second term is the fluctuation (A'fluct) contribution, 
where A^o,™ = 1/2 - ^^(k) tanh[/3S„(k)/2]/[2£;,„(k)] is 
the momentum distribution. The inclusion of A'fluct is 
very important near the critical temperature, however, 
A^o niay be sufficient at low temperatures (12llllJ|. Next, 
we discuss the T = case. 

In weak coupling (max{| Ai^oU A2,o|} ^ ^f), 
the solutions of the order parameter equation are 
max{|Ai,o|, |A2,o|} ~ Se^ exp[-2 + 7rA:o/(2A:F) - 0-] and 
min{|Ai,o|,|A2,o|} ~ Se^ exp[-2 + ^A:o/(2fcF) - '/'+], 
while the number equation leads to w ep. Here 
0± = A+ ± [\\ - 1/ det A] where A± = (An ± 
A22)/(2det A), det A = A11A22 - A12A21, and A„„i = 
VnmDm are the dimensionless interaction parameters 
with Dm — MmVkm,F/{2n^) is the density of states 
per spin at the Fermi energy. On the other hand, in 
strong coupling [fJ, < and max{|Ai,o|, |A2,o|} <C < 
eo), the solution of the order parameter equations is 
fi — — eo[7rfco/(2/ci?0+) — 1]^, while the number equa- 
tion leads to |A„,oP = (SttAA^ /M„ ) | M I / (2M™ ) . Here, 



is the fluctuation contribution to f^gauss- Expressing ^0 = fcg/(2Mi) and TV™ = N^/V. Notice that the total 



A„_o in terms of its amplitude and phase A„^o = 
\Anfi\exp{iipn) shows explicity the Josephson coupling 
energy [y22|Ai,oP + V^ii|A2,oP - 2 J| Ai,oA2,o| cos((p2 - 
</?i)]/detV of flQ. When J > 0, only the 0-phase (or in 
phase) ip2 = i-Pi solution is stable. However, when J < 0, 
only the 7r-phase (or out of phase) Lp2 ~ solution is 

stable. Thus, a phase transition occurs from the 0-phase 
to the TT-phase when the sign of J is tuned from negative 
to positive values as shown in Fig. ^p. 

From the stationary condition dSo/dA*^{q) — 0, we 
obtain the order parameter equation 



On O12 \ / Ai,o 
O21 O22 / V A2,o 



= 0, 



(3) 



density of fermions is A/" = A/i -I- A/2 = {kf^p+k2^p)/{3TT^ ). 
The familiar one-band results are recovered when J — » 
upon the use of the relation between Vnn and a„„. Next, 
we analyze the T — evolution from BCS to BEG for 
identical bands {Mi — M2 — M) with zero offset (eu — 
0). For this purpose, we set ko « 256^^ and V22 = 0.001 
in units of 5.78/ (MkpV) [or l/(fcFa22) ~ -3.38], and 
analyze two cases. 

In the first case, we solve for /i, |Ai^o| and |A2,o| as a 
function of Vu [or l/(fcFaii)], and show A„.o in Fig.l^t- 
for fixed values of J. The unitarity limit is reached at 
Vii w 1.0132^22 [or l/(A;F|aii|) w 0] while /x changes 
sign at « 1.0159 [or l/{kpaii) « 0.69]. The evolution 
of [Ai^oj is similar to the one band result where it 
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FIG. 2: Plots of (a) order parameter amplitude |A„,o| (in 
units of ep), and (b) fraction of fermions Nn/N versus Vn [in 
units of 5.78/(MA;fV)] and versus l/ikpan) for J = O.OOIV22 
(hollow squares) and J = O.OOOIV22 (solid squares). 



grows monotonically with increasing Vn. However, the 
evolution of IA2.0I is non-monotonic where it has a hump 
approximately Vn ~ 1.0155V22 [or l/{kpaii) ~ 0.58], 
and it decreases for stronger interactions until it vanishes 
(not shown). In Fig.[2l3, we show that both bands have 
similar populations for Vu ~ ¥22- However, as V11/V22 
increases, fermions from the second band are transferred 
to the first, where bound states are easily formed and 
reduce the free energy. 
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FIG. 3: Plots of (a) chemical potential /i (in units of 
ep), and (b) fraction of fermions Nn/N versus J [in units 
of 5.78/(MA:fV)] for Vn = 7V22 where 7 = l(dotted 
lines), 1.010 (solid squares), 1.014 (hollow squares), 1.016 
(crossed lines) and 1.030 (solid lines); or l/(fcFaii) ~ 
-3.38, -0.81, 0.20, 0.70 and 4.17, respectively. 

In the second case, we solve for /x, jAi^ol and |A2,o| as 
a function of J, and show /i in Fig. for fixed values 
of Vii/V22- The order parameters |Ai,o| and |A2,o| grow 
with increasing J (not shown). In Fig. we show that 
the band populations iVi and iV2 for several values of 
7. Notice the population imbalance and the presence 
of maxima (minima) in A^i (iV2) for finite J, which is 
associated with the sign change of ^ shown in Fig. 
When Vii > V22, there is population imbalance even for 
J = 0, because atom pairs can be easily transferred from 



the second band to the first until an optimal Jo is reached. 
Further increase in J produces also transfer from the first 
band to the second leading to similar populations for J S> 

Jo. 

Next, we discuss the low energy collective excita- 
tions at r = 0, which are determined by detF^^(g) = 
0. The phase-only collective excitations in weak and 
strong couplings lead to a Goldstone mode u>^(q) = 
w^jqj^ characterized by the speed of sound = 
{Divf + D2V2) / {Di + D2); and a finite frequency mode 



w^(q) 



jqj^ characterized by the finite frequency 



wl = 4a|gi2Ai,oA2,o|\/^WR(^i+^2)/(£'i£'2) and the 
speed = (i^iwl + D2vf)/{Di + 02). In weak cou- 
pling, a = 1 and w„ — w^.f/a/S, while a = l/ir and 
Vn = |A„^o|/V8M„|/i| in strong coupling. Here, Vn,F 
is the Fermi velocity. It is also illustrative to analyse 
the eigenvectors associated with these solutions in weak 
and strong couplings. In the limit of q — > 0, we ob- 
tain {61,62) oc (|Ai_oUA2,o|) for the Goldstone mode 
corresponding to an in-phase solution, while {61,62) oc 
{D2\Aifi\, — Z3i|A2,o|), for the finite frequency mode cor- 
responding to an out-of-phase solution. Our findings gen- 
eralize Leggett's weak coupling results ;15. 16]. 

Next, we discuss two band superfluidity near the crit- 
ical temperature Tc, where ]Ai.o] ~ 1^2,0] 0. For 
T = Tc, the order parameter equation reduces to 



detO = Oii022 -Oi202i =0, 



(5) 



and the saddle point number equation A^o = 
Sk n corresponds to the number of unbound 

fermions, where nF{x) = l/[exp(/3a;) -|- 1] is the Fermi 
distribution. While A^o is sufficient in weak cou- 
pling, the inclusion of iVfluct is crucial in strong cou- 
plings, and can be obtained as follows. Near T — 
Tc, the fluctuation action S'fluct reduces to S'fluct = 
(/3/2) E„,„,™ ^™U9)A;(<z)A,„(g) where 



1 - nF(6i+) - nF{S.n-) 



k 



in+ + in 



r4k)l' (6) 



corresponds to the fluctuation propagator of band n, 
L~^rnil) = 9nm, and ^„± = ^„(k ± q/2). Thus, the re- 
sulting action leads to fifiuct = X]g ln[det L"-^(g)//3^]. 
where det'L~^{q) = L^^ {q)L:^^ {q) - 312521]- No- 
tice that, detL~-'^(0) = also produces Eq. 
which is the Thouless condition. After the analytic 
continuation ivi — > w + lO^, we expand L~^{q) 
to first order in w and second order in q such 
that L-^{q) = a„ -f- c^imj I {'^M^) - dnW. 

The time-independent coefficients are given 
by a„ = -gun - Ek^n|rn(k)]V[2?«(k)] and 



c^f = Ek{^n'5y/[8a(k)] - /?r„<5„-/[166.(k)] 
/?2x„y„A;ifc,/[16M„e„(k)]}]r„(k)]2, where X„ 



tanh[/3^„(k)/2] and y„ 



sech"[/3C„(k)/2]. 



No- 



tice that, = CnSij is isotropic for the s-wave 



4 



considered here. The time-dependent coefficient 
has real and imaginary parts, and for the s-wave 
case is given by d„ = ^k^"|r„(k)|V[4C^(k)] + 

i{Pn /8)eoDn\/Jjj£FQifJ.) / {sq + fJ-), where <d{x) is the 
Heaviside function. For completeness, we present the 
asymptotic forms of a„ , c„ and dn ■ In weak coupling 
{p w ep), we find a„ = -g„„ + Dn[\n(T/Tc) + c„ = 
7eFD„C(3)/(12r27r2)^ and d„ = i?„[l/(4e;^) + z/(8T,)], 
where ({x) is the Zeta function, and Tc is the physical 
critical temperature. In strong coupling (eo ^ \n\ ^ Tc), 
we find a„ = -gn n - TTDneQ/[2^/e^{ y^\ + y^)]: 
c„ = 7rL)„/(16\/|^|eF) and d„ = vri:)„/(8-v/|At|eF). 

In order to obtain rifluct, there are two contributions, 
one from the scattering states and the other from poles 
of L(q). The pole contribution dominates in strong cou- 
pling. In this case, we evaluate deth~^{q) = and find 
the poles u;±(q) = A+ + B+|qp ± -I- B_|q|2)2 + 
512521/(^1^2)]^/^, where A± ^ (01^2 ± a2di) / {2did2) 
and B± = (M2d2Ci ± MidiC2)/(4MiM2did2)- Notice 
that, when J ^ 0, we recover the limit of uncoupled 
bands with w„(q) = fln/rfn + |qpc„/(2M„(i„). In the 
q — > limit, when J > (J < 0), the eigenvec- 
tors [A|(0), A2(0)] = [512, ai — diw±{0)] correspond to 
an in-phase (out-of-phase) mode for w+{q) and an out- 
of-phase (in-phase) mode for w_(q). Thus, we obtain 
^Ruct = (l//3)Z]±,gln[/3(wf - w{q))] which leads to 



+ gn#mA„(x) = 0, 



(8) 



fluct 



2^^j— «B[w(q)]. 



(7) 



For sufficiently strong couplings, dw±{q) /dfi = 2, and 
the poles can also be written as w±(q) = —IJ-b,± + 
|qp/(2MB_±), where /iB,± is the chemical potential 
and Mb,± is the mass of the corresponding bosons. 
In the case of identical bands with zero offset, ci = 
C2 = c and di — d2 — d, ^JLb,± = — [ai -I- 02 ± 
\/{ai - a2)2 + 4512321]/ (2d) and Mb,± = 2M in strong 
coupling. Notice that the -I- bosons always condense first 
for any J (independent of its sign) since ^ib,+ — > first. 
Next, we analyze Tc in weak and strong couplings. 

In weak coupling, solutions to the order parame- 
ter equation Eq. © are Tc,^ = (Sep/Tr) exp[7 —2-1- 
■Kko/{2kp) — (j)±], while the number equation Eq. Q 
leads to fi « ep- On the other hand, in strong cou- 
pling, the solution of the order parameter equation is 
fi = —eo[Trko/{2kp(j)^-) — 1]^, while the number equa- 
tion Af/2 = J\fB,+ +J^B.- with 7Vb,+ > Ub- leads to 

Tc,+ = T:{MB,+/m/2)iMB,+^/M^)]Y^\ since the + 
bosons condense first. Notice that, the physical critical 
temperature is Tc = max{Tc^+, Tc__} = Tc_+ for any J. 
Therefore, Tc grows continuously from an exponential de- 
pendence on interaction to a constant BEC temperature. 
Next, we obtain the TDGL equations for T « Tc, 



in the real space x — (x, t) representation. The coef- 
ficient of the nonlinear term 6„ — X^k {"''^"/[^CnC^)] ^ 
/3F„/[8^,2j(k)]} |r„(k)|'^ can also be obtained analyti- 
cally in weak 6„ = 7D„C(3)/(8Tc tt^) and strong 6„ — 
TrDn/ {■i\^J.\\/2\^\ep) coupling limits. In weak coupling, 
Eq. ||SJ| reduces to the coupled GL equations of two BCS 
type superconductors. However, in strong coupling, it 
is more illustrative to derive TDGL equations in the ro- 
tated basis of -I- and - bosons ($^,$^) = (Aj,A^)Rt, 
where R is the unitary matrix that diagonalizcs the linear 
part of the TDGL equations. In this basis, Eq. ||SJ) re- 
duces to generalized GP equations of and bosons 
showing explicitly terms coming from density-density in- 
teractions such as U±±\^±\'^^± or C/±ip|<i>± In the 
case of identical bands with zero offset, this leads to 
U++ = U— = h/{2d^) and [/+„ = = 36/(2^^) 

as the repulsive density-density interactions. 

In conclusion, we showed that a quantum phase tran- 
sition occurs from a 0-phase to a 7r-phase state depend- 
ing on the relative phase of the two order parameters, 
when the interband interaction J is tuned from negative 
to positive values. We found that population imbalances 
between the two bands can be created by tuning intra- 
band or interband interactions. In addition, we described 
the evolution of two undamped low energy collective ex- 
citations corresponding to in-phase phonon (or Gold- 
stone) and out-of-phase exciton (finite frequency) modes. 
Near the critical temperature, we derived the coupled 
Ginzburg-Landau (GL) equations for a two-band super- 
fiuid, and showed that they reduce to coupled Gross- 
Pitaevskii (GP) equations for two types of composite 
bosons in strong coupling. 
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[a„ + 6„|A„(a;)|^ 
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